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SUMMARY 

This report addresses techniques for modeling complex, unclasslcal 
contact-friction problems arising In solid and structural mechanics. A 
constitutive modeling concept Is employed whereby analytic relations between 
Increments of contact surface stress ( 1 . e . , traction) and contact surface 
deformation (l.e., relative displacement) are developed. Because of the 
Incremental form of these relations, they are valid for arbitrary load- 
deformation histories. The motivation for the development of such a constitu- 
tive law Is that more realistic friction Idealizations can be Implemented In 
finite element analysis software In a consistent, straightforward manner. Of 
particular Interest In this report Is modeling of two-body (l.e., unlubricated) 
metal-metal, ceramic-ceramic, and metal-ceramic contact. Interfaces Involving 
ceramics are of engineering Importance and are being considered for advanced 
turbine engines In which higher temperature materials offer potential for 
higher engine fuel efficiency. 


1. INTRODUCTION 

Material Interfaces are common In mechanical systems and often have a 
substantial Influence on structural response. The behavior of a material 
discontinuity Is complex and Involves frictional sliding, possible contact 
surface separation and various types of surface damage. Because of the non- 
linearity of such behavior, numerical methods such as the finite element method 
(FEM) are popular for the analysis of problems Involving discontinuities. How- 
ever, solution methodologies are still not advanced to the point where contact- 
friction capabilities are Included In general purpose FEM programs and In most 
analyses, special purpose analysis programs are used. 

Because of the difficulty In analyzing contact-friction problems and our 
rather poor quantitative understanding of frictional phenomena, virtually all 
of the analyses reported In the literature have employed a very simple and 
Idealized form of friction law which can be attributed to Coulomb's original 
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work published In 1781. In Its "modern" form, this law postulates that the 
magnitude of the tangential component of the surface traction, or simply the 
tangential stress, will not exceed the product of the normal component of the 
traction, or normal stress, and a constant called the coefficient of friction. 
For states of tangential stress below the friction limit, various types of 
behavior such as stlck-lmpenatrablllty or penalized deformablllty are common. 
Remarkably, for many problems, Coulomb's law Is accurate and Is sufficiently 
faithful to reality for engineering analysis. In fact, much of the research 
In tribology since Coulomb's work has been directed at obtaining a qualitative 
understanding of why and under what circumstances this law Is accurate. How- 
ever, It has been established by experimentation that true friction behavior 
Is considerably more complicated than Coulomb's Idealization, although for 
engineering analyses, It Is still an open question when such detail Is 
warranted. 

In this paper, we employ a constitutive modeling concept to develop 
analytic relations between Increments of contact surface stress (l.e., trac- 
tion) and deformation (l.e., relative surface displacement). The development 
begins by distinguishing between macrostructural and mlcrostructural features 
of a contact surface. Through macrostructural considerations and the assump- 
tion that contact surface deformations consist of reversible (elastic) plus 
Irreversible (plastic) parts, the general form of the constitutive law Is 
obtained. Mlcrostructural considerations permit the specialization of the 
constitutive law for various applications by allowing the Inclusion of special 
frictional phenomena such as adhesion, mlcroslldlng effects on coefficient of 
friction and surface damage. In particular, the behavior of unlubricated 
metal-metal, ceramic-ceramic, and metal-ceramic contacts will be discussed. 

Because of the simplicity of Coulomb's law, most analysis programs employ 
ad hoc Implementations of the contact-friction conditions which usually are not 
amenable to Implementation of more general friction laws. Our objective Is the 
development of a general, explicit. Incremental constitutive law that leads to 
straightforward Implementation In finite element software. In Section 4 of 
this report, a simple finite element and solution procedure for plane contact 
problems Is discussed. 


NOMENCLATURE 


A actual contact area 

A q available or macroscopic contact area 

B j denotes body number I where 1=1,2 

E^ Interface stiffnesses where 1 = t, n and j = t, n 
e^ unit vector In coordinate direction 1= t, n 

F slip function 

G slip potential 

g^ relative discontinuity displacement In coordinate direction 1 = t, n 
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H hardening or softening parameter 

L length of contact finite element 

Nj . finite element shape function for node I = 1, 2, 3, 4 

n surface Index for finite element n = (-1 surface 1-2 

1+1 surface 3-4 

p arbitrary contact point between two bodies 

sgn(ot) denotes the sign of sgn(<*t) = /+! If a t > 0 

l-l If ot < 0 

t variable denoting tangential distance along finite element 

u T displacement vector of a point associated with body B T ; I * 1,2 

~l l 

W p plastic work 

asperity angle 

at time step for dynamic analyses 

n exponent used In friction law given by equation (3.7) 

V scalar giving magnitude of plastic slip 

y friction coefficient 

\ ratio of actual contact area to available contact area 

Interface stress; 1 = t, n 
Matrix Symbols 

matrix relating relative displacements to finite element node 
displacements 

vector of nodal forces 

element matrix Including nonllnearltles 

shape function matrices 

transformation matrices 

vector of nodal displacements 

Subscripts 

1 ,j denotes coordinate direction t or n 

I denotes body number B-| and 83 or finite element node numbers 1 

through 4 
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B 

t 

k 6P 
M, N 

X, 


K denotes sliding direction to right or to left 

n,t denotes directions normal and tangential, respectively, to the 

macroscopic contact surfaces 

y used In equation (3.8) denoting contact yield stress with a certain 

amount of constrained plastic deformation 

denotes a vector or matrix 

Superscripts 

e,p denotes elastic and plastic, respectively 

j denotes state of stress and deformation at Iteration j for 

Incremental static analysis or time jAt for dynamic analysis 

(•) denotes the time rate of change 


2. CONSTITUTIVE LAW - MACROSTRUCTURAL CONSIDERATIONS 

Shown In figure 1 Is a macroscoplcally smooth contact surface with local 
tangential and normal coordinate directions t and n, respectively, with 
origin at point p~ which Is affixed to body B-| . Any roughness that may 
be present on the contact surface Is not shown In figure 1 and will be consid- 
ered as mlcrostructural features which are discussed In Section 3 of this 
paper. The only requirement that the macroscopic contact surface must possess 
Is that It be sufficiently smooth so that the surface tangent and normal are 
not Ill-defined. 

The kinematic variables to be used In the constitutive law are the rela- 
tive surface displacements In the tangential and normal directions which are 
defined as 


91 = (Ai2 - “I) * ®1 (2*1 ) 

where ui Is the displacement of the point p associated with body I and 
e^ Is a unit vector In coordinate direction 1 = t,n. The stresses (In 
proper terminology, these are tractions, but we adopt the more conventional 
nomenclature of stresses) the Interface supports at point p are denoted by 
and the convention that compressive stresses are negative Is used. Our 
objective Is the development of an explicit relation between Increments of 
stress and deformation. Because the relation Is Incremental, It Is valid for 
arbitrary load and deformation histories that can Involve unloading and subse- 
quent reloading ( 1 . e . , changes In the direction of frictional sliding). 

A basic assumption Is that the deformation can be addltlvely decomposed 

Into 


g 1 - g* + 9? (2.2) 

where superscripts e and p denote the elastic (recoverable) and plastic 
(Irrecoverable) parts of the deformation. There exists experimental evidence 
for almost every class of friction problem that has been carefully studied 
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supporting this decomposition. Furthermore, It has the added advantage of 
leading to a more convenient numerical Implementation compared to friction 
Idealizations In which a "stick" condition precedes frictional sliding. 

Because the stress that an Interface supports must be reversible (l.e., 
stresses must go to zero upon unloading). It can be related to the elastic part 
of equation (2.2) only. The simplest relation possible Is the linear law 

V E 1j9j < 2 - 3 > 

where the are constant Interface stiffnesses, superposed dots denote 

time differentiation and the summation convention Is applied to repeated 
Indices. Although nonlinear E^j are possible. It Is not known If such 
detail Is warranted. Based upon physical considerations. It Is appropriate to 
take E^n = E n t = 0 so that changes of stress In the tangential and normal 
directions are unrelated to changes of deformation In the normal and tangential 
directions, respectively. For example, expanding the second of equations (2.3) 
gives o n = E n ^gt + E nn g n . If changes In tangential elastic displacement 
are not to give rise to changes In normal stress, then E n t =0- A similar 
argument can be made showing Etn 3 0. 

It Is necessary to postulate a rule for the plastic deformation. By 
assuming that 

(1) A slip potential (scalar valued) function, F, can be defined such that 
negative F corresponds to nonslldlng states of stress, zero F corresponds 
to sliding and positive F Is not possible, and 

(2) Increments of stress are linearly related to Increments of 
deformation. 

It can be shown (see ref. 1) that a slip rule, which Is analogous to the flow 
rule employed In continuum plasticity. Is the appropriate law for the plastic 
deformations 


0 If F<0 or F < 0 (unloading) 

. 3G , (2 * 4) 

X ~ If F = f = 0 (loading) 

where, G Is the slip potential whose gradient gives the direction of the slip 
and x Is a nonnegative scalar that gives the magnitude of the slip. When 
F = G, the friction Is "associated" and when F * G the friction Is "non- 
assoclated". Only under unusual circumstances Is friction associated. 
Unfortunately, the more typical case of nonassoclated friction results In a 
nonsymmetrlc tangent material matrix which poses difficulties for many computer 
solution schemes; methods of reducing these difficulties will be mentioned In 
Section 4 of this report. 

Combining equations (2.1) through (2.4) provides the general form of the 
contact-friction constitutive law 
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(2.5) 


• c e P A 

°1 = E U 9 J 

where for purely elastic response or unloading (l.e., F < 0 or f < 0), 
respectively). 


E ep - E 

E 1j " 


and for Imminent slip (l.e., F = f = 0) 


(2.6(a)) 


Eff.E 


it- t P ifi. 
3 a 1q pi 3 a 

p q 


1j ~ u 1j a£_ E 36 H 
ao p pq a 0q 

where H Is a hardening or softening parameter given by 


(2.6(b)) 


H . a. ~- 

* 8 °l 

In which W p = oi g£ Is the plastic work. If no hardening or soften- 
ing behavior Is permitted, then H = 0. 


(2.6(c)) 


Remarks 

1. Completion of the constitutive law requires the specification of the 
slip function F and slip potential 6. These are determined by considering 
the microstructure of the material Interface and will be considered In Section 
3 of this report. 

2. Equation (2.5) Is an explicit. Incremental relation that Is In a form 
convenient for numerical Implementation and the analysis of problems with 
arbitrary load and deformation histories. 

3. This constitutive law enforces compatibility (l.e. Impenetrability 
constraint preventing contacting bodies from penetrating one another) through 
the use of a “penalty' 1 stress that Is proportional to the violation of compat- 
ibility. Thus, when bodies are In contact, gfi < 0 (bodies are allowed to 
Interpenetrate slightly) and the penalty stress Is o n = E n n9fi* 


3. CONSTITUTIVE LAW - MICROSTRUCTURAL CONSIDERATIONS 

The general framework of the constitutive law given by equations (2.5) and 
(2.6) Is applicable to a variety of contact problems. These varied applica- 
tions are obtained by formulating the slip function F and slip potential 6 
such that Important physical features of the contact problem are accounted for. 
In this section we demonstrate several examples. 
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The first example to be considered Is Illustrative, simple, and Identical 
to the penalized form of the classic Coulomb friction law (ref. 2). We ideal- 
ize the material discontinuity as smooth with no roughness whatsoever. 
Coulomb's law of friction requires 


|o t ! - ~ va n 

The slip function and slip potential for this case are 

F = I + v<» n 


(3.1) 


(3.2) 

These surfaces are shown In figure 2. Equations (2.5) and (2.6(b,c)) become 


G = k t l 


a t = -(sgn(o t ) v E n ) g p 


(3.3) 


^n^n 


(3.4) 


where sgn(c^) denotes the sign of and y Is the coefficient of 
friction. Equation (3.4) Is seen to be the correct relation between rates of 
normal stress and normal relative displacement and equation (3.3) provides the 
relation between the tangential stress and normal deformation rates where 
E n g n Is recognized as the rate of change of normal stress and sgn(a^) 

Insures that the tangential stress rate has the appropriate sign. 

Rough, dllatant material discontinuities . - Prototypical material Inter- 
faces are microscopically very rough; two extreme surface profiles are shown 
In figure 3. The situation depicted In figure 3(a) has a rough surface with 
close Initial mating and Is characteristic of Interfaces, or cracks, that 
propagate through an Initially continuous medium. An Important feature of 
these types of discontinuities Is dllatancy: contact surface separation during 

relative tangential motion due to the asperity surfaces of one body riding up 
on those of the other. Such a discontinuity Is typical of crack surfaces In 
polycrystalline and aggregate materials and the effects of dllatancy may be 
Important. For example, consider a plane tortuous discontinuity with a sharp 
crack front In a material that Is loaded In plane shear or antiplane shear. 

If the crack was not tortuous then these loading conditions would give rise to 
mode II and III crack growth, respectively. However, because of the rough, 
closely mated discontinuity, mixed mode crack growth (addition of mode I) may 
result due to dllatancy. A simplified model for the surface shown In 
figure 3(a) was developed In references 3 and 4 with the following slip func- 
tion and potential 

F = |o Sin + a. COS a„| ♦ y(o COS a„ - a. sin a„) (3.5) 

n k t k n k t k 

G = |c sin a + a. COS a I (3.6) 

n In t In 

where a* Is an average asperity orientation with respect to the macro- 
scopic Interface tangent and K Is either R or L to denote if sliding Is 
In the right hand or left hand direction, respectively. Equations (3.5) and 
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(3.6) were obtained by assuming that simple Coulomb friction as given by 
equations (3.3) and (3.4) Is active on each asperity surface. The surface 
F = 0 Is shown In figure 4; the step discontinuity In F occurs when the 
asperities have overridden one another and the discontinuity behavior becomes 
nondllatant and Is assumed to be governed by the simple Coulomb friction law 
given by equations (3.3) and (3.4). More realistic friction laws, such as 
equation (3.7) discussed In the following subsection, are possible for dllatant 
surfaces, however, this might correspond to the modeling of second order 
effects since for rough, closely mated contact problems, dllatancy Is probably 
responsible for the behavior that Is observed. 

Rough, nondllatant material discontinuities . - The situation depicted In 
figure 3(b) Is a rough Interface with unmatched surfaces so that contact takes 
place at only a small number of asperity peaks. Such a discontinuity Is 
characteristic of man-prepared surfaces that come Into contact. For engineer- 
ing purposes, these Interfaces are nondllatlng with sometimes complex fric- 
tional behavior. The friction of such contacting solids Is extremely sensitive 
to a number of factors Including environment, roughness, presence of contami- 
nating film layers, sliding history, and temperature. Thus, It Is very diffi- 
cult to make generalities about frictional behavior, although In this paper we 
will attempt to do this with the understanding that there are numerous excep- 
tions. In postulating models for this behavior, we draw heavily upon the work 
of Courtney-Pratt and Eisner (ref. 5), Bowden and Tabor (ref. 6), Kragelskll 
(ref. 7), and Buckley (ref. 8). Experimental methods for the determination of 
necessary friction law parameters are not discussed and may pose some open 
questions requiring additional research. 

In many respects, metals, glasses and ceramics have similar friction and 
wear behavior. In this paper, we will attribute two-body (unlubricated) fric- 
tion to three primary sources: (1) plasticity of the surface film layer, (2) 

adhesion due to the molecular attraction of bodies In very close proximity, and 
(3) plowing due to the penetration of prominent asperities through a surface 
film layer. If any. Into the adjacent material. Usually only one of these 
sources of friction will be active at any given time although sometimes sources 
(1) and (3) or (2) and (3) can exist simultaneously. To represent this behav- 
ior the following friction law Is considered 


(3.7) 


The first term of equation (3.7) represents the basic frictional resis- 
tance of the film (contaminant) layer to plastic deformations. Surface films 
are almost always present In noninert environments such as air. The frictional 
properties and wear resistance of most materials Is extremely sensitive to 
environment and with the notable exception of glass, most materials exhibit 
reduced friction and wear In film-producing environments compared to their 
behavior In a vacuum In which surface films are not present or are rapidly 
destroyed without regeneration. Compared to their frictional resistance In a 
film-producing environment, ceramics and glasses show up to a 30 percent 
Increase In a vacuum, while the Increase for metals Is significantly greater 
at up to two orders of magnitude (complete seizure). The reason for this 
marked behavior is a change In the source of friction. In the presence of a 
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surface film, the second and often times third terms of equation (3.7) are 
negligible while the opposite Is true In an Inert environment (ref. 8). 


The second term of equation (3.7) represents frictional resistance due to 
molecular attraction between bodies In very close contact. Thus, the larger 
the number of asperities In Intimate contact, l.e., the larger the actual area 
of contact, the greater the adhesion. In equation (3.7), 5 Is the ratio of 
the actual contact area. A, to the available, or macroscopic area, A 0 . To 
quantify this, we consider the model of asperity contact shown In figure 5 In 
which the state of stress at the asperity Is assumed to be sufficiently high 
to cause local yielding. Assuming a maximum shear stress yield criterion 
provides 


5 


2 




(3.8) 


where ay Is the yield stress, although It must be Interpreted In a slightly 
different sense than Its usual definition as the uniaxial yield stress since 
behavior of asperities Involves local large deformation yielding. Rather, It 
can be Interpreted as the stress necessary to produce a certain amount of 
constrained plastic deformation which Is significantly higher than the Initial 
yield stress (ref. 9, pg. 335). This constant and vattr are determined 
from experiments In which conditions are selected to eliminate the other 
sources of friction appearing In equation (3.7). 


Where surfaces are atomically clean, this source of friction can be very 
substantial and particularly with metals, the frictional resistance from 
adhesion and plowing can approach full seizure. Anything that tends to disrupt 
atomically clean surfaces such as surface films or contaminants resulting from 
an air environment decreases v a ttr an£ * ma k es the affects of vfiim more 
significant. 

The last term of equation (3.7) represents friction due to the plowing of 
asperities of one body through those of the other. This source of friction 
becomes more prominent at higher compressive loads but Is diminished by the 
presence of film layers which add a zone of softer material that must be 
penetrated. For ceramics, frictional resistance Is generally load Independent 
except when a transition Is made when asperities penetrate a surface layer and 
begin plowing. The normal stress at which seizure occurs, o S e1ze» alon 9 with 
n and ypiow must & e determined from experiment. This source of friction 
Is prominent when contact Involves materials with disparate hardnesses. For 
example, ceramic-ceramic and glass-glass contact Is controlled by film layer 
and adhesive friction with no plowing. However, when metal contacts ceramic, 
the metal Is much softer and the ceramic will plow offering additional 
resistance. 


Although the constitutive law presented In this paper does not model wear, 
It Is an Important phenomenon that can Influence frictional behavior. Because 
of adhesion, there often Is a transfer of one material to another thus affect- 
ing a change In surface profile and characteristics. This type of wear Is most 
pronounced for metals In contact with ceramic and glass (ref. 8). With metal- 
ceramic contact systems, metal transfers to ceramic and eventually the contact 
behaves as a metal-metal system. With metal-glass In air, metal also transfers 
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to glass although when In a vacuum, the opposite Is generally the case; the 
reason for this anomalous behavior Is that the fracturing resistance of glass 
Is strongly water moisture dependent. 

Anisotropy of frictional sliding and wear has been experimentally observed 
for many single crystal metal and ceramic contact systems. For this case, a 
version of equation (3.7) can be employed for each of the two principal tangen- 
tial directions (directions displaying maximum and minimum values of friction). 
Slip functions and slip potentials for anisotropic friction are presented In 
references 10 and 11 . 


4. FINITE ELEMENT IMPLEMENTATION 

In this section we review the development of a simple two-dimensional 
linear displacement finite element for contact-friction problems (ref. 2). The 
element Is suitable for modeling contact between regions that are discretized 
by, for example, constant strain triangles or bilinear displacement quadri- 
laterals. The macroscopic contact surface shown In figure 1(b) Is discretized 
Into elements. One such element Is shown In figure 6 with a natural (t,n) 
coordinate system having origin at the center. The element has zero thickness 
In the undeformed state but Is shown separated for clarity. 

Kinematics . - The element has four nodes and each node has two degrees of 
freedom corresponding to displacements ut and u n , tangent and normal to 
the plane of the surface, respectively. The length of the element Is L. The 
displacement of any point on surface 1-2 or surface 3-4 can be expressed In 
terms of nodal displacements by 


u 1 ( t,n) = N i u 1I I = 1 ^2,3,4 (4,1) 

where a lower case subscript denotes a vector component, upper case subscript 
I denotes a node number, and the shape functions are 

"i ■ H 1 - r) d - "> 

"2 ■ 4 l 1 ♦ F) <’ - "> 

(4.2) 

N 4 ■ H 1 - f) (1 ♦ "> 

where -L/2 < t < + L/2, n = -1 for surface 1-2 and n = +1 for surface 3-4; 
values of n such that -1 < n < +1 have no meaning since the Interface Is 
assumed to have Infinitesimal thickness when In contact. In matrix notation, 
equation (4.1) becomes 

u(t,n) = Nu (4.3) 
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where 


u(t,n) = [u t (t,n), u n (t,n)] 


u * [u tl , u t2 , u t3 , u t4 , u nl , u n2 , u n3 , u n4 ] 


( 4 . 4 ) 


N = 


N 


1 


N 2 N 3 N 4 


0 

N, 


'1 "2 N 3 ” 4 J 

Superscript T denotes vector or matrix transposition. 

The relative tangent and normal contact displacements are given by 

g^(t) = u^t, +1) - u^t, -1) 


0 

N. 


( 4 . 5 ) 


The distribution of the relative displacements, g^ ( t ) , can be related to the 
node relative displacements by 


a(t) = Hg 


( 4 . 6 ) 


where 


a T (t) = [9 t (t), g n (t)] 


a - [g t r g t2’ g m’ V 1 


m = 


M 1 M 2 0 0 


0 0 M 1 M 2 


The shape functions M are related to the shape functions N by 


( 4 . 7 ) 


H = NT 1 


( 4 . 8 ) 


where 


~1 


1 

0 

0 

1 

0 

0 

0 

.0 


0 

1 

1 

0 

0 

0 

0 

0 


0 

0 

0 

0 

1 

0 

0 

1 


0 

0 

0 

rt 

u 

0 

1 

1 

0_ 


( 4 . 9 ) 
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Furthermore, the relative node displacements are related to the total node 
displacements by 


S = I2li 


(4.10) 


where 


I 2 



0 

-1 

0 

0 


0 

1 

0 

0 


10 0 0 
0 0 0 0 
0-100 
0 0-11 


0 

0 

1 

0 


(4.11) 


Combining equations (4.6), (4.8), and (4.10), the distribution of the 
relative displacements are related to the total node displacements by 


g(t) = Bu 


(4.12) 


where 


S - HXi I 2 


(4.13(a)) 


If the matrix product In equation (4.13(a)) Is expanded, we obtain 
(1-t) -(Ht) (1+t) (1-t) 0 0 

0 000 -(1-t) -(1+t) (1 

(4.13(b)) 


•A 


0 0 1 
l+t) (l-t)J 


EquI librium . - To equilibrate the distribution of the contact stresses, 
o(t), with the slave nodal (concentrated) forces, £, we will require the work 
of the distributed stresses to be equal to the work of the nodal forces. Thus, 


A/2 

J-L/2 


d£ T (t) t 2 J (t) + do(t) ] dt = du T f J+1 


(4.14) 


where 


~ = t f tl* f t2’ f t3’ f t4’ f nl * f n2’ f n3’ f n4 J 


(4.15) 


and superscript j denotes, In the case of a dynamic analysis, time jAt 
where At Is the time step and In the case of an Incremental static analysis, 
j denotes the state of stress and deformation at Iteration number j. 


In 


Combining (4.14) with (4.12) and noting the arbitrariness of du results 




1 C U2 T 
= f J + I B* do 

J -U2 ~ 


do(t) dt 


(4.16) 
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where 


= f B T c J (t) dt (4.17) 

•U/2 

or In other words, the Increment of d£; d£ = - fJ, Is related to the 

Increment of the stresses by 


f U2 

df = B 1 do(t) dt (4.18) 

•U/2 

where B and da are given by equations (4.13) and (2.5), respectively. 

Solution Schemes . - For transient finite element analysis, the Implementa- 
tion of the constitutive law presented In this paper Is particularly straight- 
forward If the time Integration scheme Is explicit. Two popular explicit 
Integrators are the central difference and the explicit Newmark methods 
(ref. 13). The advantage of explicit methods Is that at each time step the new 
displacement configuration Is determined completely In terms of historical 
Information. This enables equation (4.18) to be evaluated using stresses based 
upon historical displacements. Furthermore, If a lumped (diagonal) mass matrix 
Is employed, time stepping can be accomplished without solving systems of non- 
linear simultaneous equations. Explicit methods are easy to program and show 
good convergence for nonlinear problems but are conditionally stable and often 
permit only small time step sizes. Nevertheless, for short and moderate dura- 
tion transients, these methods can be very good. Implicit methods, such as 
trapezoidal rule and the Implicit Newmark family, compute a new displacement 
configuration using historical Information and derivatives of the new (unknown) 
displacement. Thus, these methods require the solution of a nonlinear system 
of simultaneous equations at each step which can entail considerable diffi- 
culty. The attractive feature of these methods Is that they have good stabil- 
ity properties that permit the use of time steps that are larger than those for 
explicit methods. However, these methods are more difficult to program and 
convergence for large time steps with the constitutive law presented In this 
paper Is an open question. 

For static analysis, an Incremental solution scheme can be obtained by 
combining equation (4.18) with equations (4.12) and (2.5) to obtain 

df = k ep du (4.19) 

where the Instantaneous (tangent) stiffness matrix Is 

L/2 

k ep = f B T E ep B dt (4.20) 

•U/2 

As previously mentioned, £ e P and hence, k®P Is unsymmetrlc for 
frictional response. In addltTon, It Is possible that It may occasionally be 
singular. Thus, approximate tangent stiffness matrix methods (refs. 14 and 15) 
( 1 . e . , approximate Newton-Raphson methods) that use symmetric matrices may be 
more robust and economical. 


13 


SUMMARY AND CONCLUSIONS 


A constitutive model for complex, unclasslcal contact-friction problems 
has been presented for applications to dllatant and nondllatant surfaces. 
Surface deformation In the normal and tangential directions are decomposed Into 
elastic (recoverable) and plastic (Irrecoverable) components. Resultant normal 
and tangential stresses are determined from the deformations through a general 
elasto-plastlc material matrix that accounts for several Important physical 
features of the contact problem, namely surface roughness, relative mating of 
contact surfaces, and an appropriate surface friction law. For the unlubri- 
cated two body contact problem considered here, a three term general friction 
law was formulated Including plasticity of a contact surface film layer, adhe- 
sion due to molecular attraction, and plowing due to penetration of contacting 
surface asperities. Application of this law to metal-metal, ceramic-ceramic, 
and metal-ceramic contact systems was discussed. Under general conditions, 
usually only one or two terms of this model would actually be required to 
adequately describe the contact. 

Finally, a linear displacement contact Interface finite element that Is 
compatible with constant strain triangle or bilinear displacement quadrilateral 
finite elements was formulated. Implementation of this element Into computer 
codes combined with the solution schemes discussed should afford sufficient 
convergence rates and solution stability. Based on this study, the following 
specific results were obtained: 

1. The constitutive law presented In this paper results In an elasto- 
plastlc material matrix for general frictional behavior and leads to straight- 
forward finite element Implementation and solution using existing 
methodologies. 

2. The Incremental form of the constitutive relation allows modelling of 
complex load-deformation histories and can account for reversals of sliding 
direction at any time. 
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(a) (b) 

Figure 1.- (a) Two body contact in 2-dimensions, (b) macrostructural contact surface 

WITH COORDINATE SYSTEM AND SURFACES SHOWN SEPARATED FOR CLARITY. 
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Figure 2. - (3) Slip surface F = 0. Slip surface and slip 

POTENTIAL FOR NONDILATANT COULOMB FRICTION. THE SURFACE 
G WHOSE NORMAL GIVES THE DIRECTION OF SLIDING, g? . IS 
ALSO SHOWN, (b) POSSIBLE STRESS-DEFORMATION RESPONSE SHOW- 
ING UNLOADING AND SLIDING IN THE OPPOSITE DIRECTION. 



(b) 


Figure 3.- Examples of rough surface profiles; (a) dilatant 

SURFACE PROFILE WITH VERY CLOSE MATING, (b) NONDILATANT 
SURFACE PROFILE WITH CONTACT AT ISOLATED ASPERITY PEAKS 
ONLY. 





Figure 4.- Slip function F = 0 for dilatant Coulomb friction. 


Figure 5.- Idealization of plastic 

DEFORMATION BETWEEN ASPERITIES IN 
CONTACT. 
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Figure 6.- 2-Dimensional linear 
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